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We study quantum correlations beyond entanglement in two-mode Gaussian states of continuous 
variable systems, by means of the measurement-induced disturbance (MID) and its ameliorated 
version (AMID). In analogy with the recent studies of the Gaussian quantum discord, we define a 
Gaussian AMID by constraining the optimization to all bi-local Gaussian positive operator valued 
measurements. We solve the optimization explicitly for relevant families of states, including squeezed 
thermal states. Remarkably, we find that there is a finite subset of two-mode Gaussian states, 
comprising pure states, where non-Gaussian measurements such as photon counting are globally 
optimal for the AMID and realize a strictly smaller state disturbance compared to the best Gaussian 
measurements. However, for the majority of two-mode Gaussian states the unoptimized MID 
provides a loose overestimation of the actual content of quantum correlations, as evidenced by its 
comparison with Gaussian discord. This feature displays strong similarity with the case of two 
qubits. Upper and lower bounds for the Gaussian AMID at fixed Gaussian discord are identified. 
We further present a comparison between Gaussian AMID and Gaussian entanglement of formation, 
and classify families of two-mode states in terms of their Gaussian AMID, Gaussian discord, and 
Gaussian entanglement of formation. Our findings provide a further confirmation of the genuinely 
quantum nature of general Gaussian states, yet they reveal that non-Gaussian measurements can 
play a crucial role for the optimized extraction and potential exploitation of classical and nonclassical 
correlations in Gaussian states. 

PACS numbers: 03.67.-a, 03.65.Ta, 42.50.Dv 
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I. INTRODUCTION 

One of the seminal findings triggering the development 
of quantum information theory is that there exist nonlo- 
cal correlations among subparts of quantum systems that 
do not emerge in a classical scenario. These nonclassi- 
cal correlations, commonly identified with entanglement, 
can be exploited to manipulate and transmit information 
in novel and enhanced ways [ij going beyond the possi- 
bilities of classical physics. Consequently, an increasing 
interest in their study has risen in recent years 0. 

Interestingly, signatures of correlations having no clas- 
sical counterpart can be traced even in separable (non- 
entangled) states, but their nature is rather different from 
entanglement [1, In fact, while entanglement is a 
consequence of the superposition principle, more general 
forms of nonclassical correlations arise essentially from 
the noncommutativity of quantum observables. When 
speaking about composite systems, separable states are 
often perceived as essentially classical. However, truly 
classical states, i.e. states which contain only classical 
correlations, represent just a subset of separable states. 
Moreover, it is also possible to show that almost all sep- 
arable states possess a finite amount of nonclassical cor- 
relations m. This has fueled a still unsettled debate, and 
an active stream of research, to decide whether separa- 
ble states containing nonclassical correlations can be also 
directly useful for quantum information tasks [B-dl ■ 

To quantitatively assess various aspects of entangle- 



ment, several entanglement measures have been adopted 
and extensively studied . In a similar fashion, more re- 
cently a zoology of indicators of nonclassical correlations 
(in separable or entangled states) have been introduced 
d, 13,11, [Tol - [l3j , among which the most popular being 
quantum discord A nice facet of discord is that it 

has an immediate characterization in information theory, 
it is endowed with operational interpretations p^ . and 
its evaluation for bipartite states of simple quantum sys- 
tems such as two-level systems (qubits) is conceptually 
straightforward, albeit technically hard. 

Studies on nonclassicality indicators are not restricted 
to finite dimensions, but have also been extended to 
quantum systems with infinite-dimensional state space, 
where correlations are encoded between continuous vari- 
ables (CVs), that is, variables with continuous spectra. 
They are represented by, e.g., modes of the electromag- 
netic field described by quadrature amplitudes. For these 
systems, a privileged role is played by the states pos- 
sessing a Gaussian Wigner function, the so called Gaus- 
sian states, as they are easy to handle both theoretically 
and experimentally [3]. Although Gaussian states 
are sometimes flagged as the "most classical" class of 
continuous-variable (CV) states because of the positiv- 
ity of their Wigner function, recent results of the non- 
classicality indicators show the opposite. Namely, a re- 
cently derived analytical form of quantum discord for 
two-mode Gaussian states [l3, 13 reveals that, contrar- 
ily to the above naive categorization, all non-product bi- 
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partite Gaussian states have a nonzero discord and so 
exhibit nonclassical correlations. 

Another nonclassicahty indicator frequently employed 
in the literature is the measurement-induced disturbance 
or MID, introduced in [l^ . A good property of MID is 
the intrinsic symmetry under swapping of the subsystems 
(unlike discord) but the main flaw is that it does not 
incorporate any optimization over local measurements, 
therefore usually returning an an overestimation of the 
actual amount of nonclassical correlations. The neces- 
sity of a more faithful nonclassicahty quantifier moti- 
vated the introduction of a new indicator called AMID 
(ameliorated measurement-induced disturbance) which is 
an ameliorated synthesis of the discord and MID, being 
both symmetric as MID and optimized as discord. It has 
been applied to characterize nonclassical correlations of 
arbitrary two-qubit mixed states 12|. The AMID is de- 
fined as the difference between the total (quantum) and 
the classical mutual information [ll|, [l^, [H, [13] , where 
the latter quantity represents the maximum classical cor- 
relations that can be extracted from a bipartite state via 
local measurements on subsystems. Here an optimiza- 
tion over all possible bi-local measurements is involved, 
which represents the nontrivial part in calculation of the 
AMID. 

In the present work, we continue the programme of 
characterizing quantum versus classical correlations be- 
yond entanglement in CV systems. Specifically, we ex- 
tend the definitions of MID and AMID to the paradig- 
matic class of two-mode Gaussian states. Contrary to 
the case of discord, which has been computed for Gaus- 
sian states restricting the optimization to Gaussian mea- 
surements only [TtI . [18| , for the nonclassicahty indicators 
considered here non-Gaussian measurements play an im- 
portant role. The MID is in fact obtained as the state 
disturbance induced by local projections onto Fock states 
(photon counting), a clearly non-Gaussian measurement. 
This can be done analytically for simple cases (such as 
pure two-mode squeezed states) and requires numeri- 
cal evaluation in more general Gaussian states. For the 
AMID, a competition between photon counting and op- 
timized Gaussian measurements takes place in the maxi- 
mization of the classical mutual information. We define a 
Gaussian version of AMID (Gaussian AMID), restricted 
to Gaussian measurements in analogy with the Gaussian 
discord [13, [11] , and provide the analytical framework for 
its computation, obtaining simple closed forms for rele- 
vant families of states, including pure states and squeezed 
thermal states. We compare this measure with MID and 
Gaussian discord, elucidating similarities and differences 
with the two-qubit case [12]. Similarly, we find that also 
for Gaussian states, in most cases, MID significantly over- 
estimates quantum correlations, and we provide evidence 
for states with nearly vanishing discord but arbitrarily 
large MID. Yet, quite interestingly, we find a finite re- 
gion of two-mode Gaussian states where MID is strictly 
smaller than the optimal Gaussian AMID, meaning that 
in those instances non-Gaussian measurements are neces- 



sarily optimal for the calculation of the classical mutual 
information and for the AMID. This subset of states in- 
cludes, surprisingly, pure two-mode squeezed states, for 
which local photon countings are found to be therefore 
"less disturbing" than local homodyne detections (which 
realize the optimal Gaussian measurements in this case), 
and allow one to extract strictly more strongly corre- 
lated classical random variables from the pure Gaussian 
states. The gap between photon counting and homodyne 
detection in the degrees of classical mutual information 
and AMID persists even in the limit of infinite squeezing. 
Such a finding is quite unexpected and certainly deserves 
further investigation with an eye on potential practical 
applications. Yet, this agrees in spirit with with a series 
of somehow related results: For example, the best oper- 
ation to clone Gaussian coherent states, as well as the 
best partial measurement achieving the optimal infor- 
mation/disturbance tradeoff for the same type of state, 
and the best measurement localizing maximum bipar- 
tite entanglement for Gaussian states, are all known to 
be non-Gaussian [2l| - [23| . Recall also that non-Gaussian 
measurements are necessary for universal CV quantum 
computation with Gaussian cluster states 24| . 



The paper is organized as follows. In Section [IT] we set 
up the notation and recall the main concepts of Gaus- 
sian states and Gaussian measurements, followed by an 
overview of the nonclassicahty measures landscape. The 
derivation of a manageable expression of MID (associated 
to local Fock projections) for two-mode Gaussian states, 
in closed form for pure states, is presented in Section [IIII 
In Section IIVI we face the evaluation of the AMID for 
arbitrary two-mode Gaussian states. We bound AMID 
from above by the minimum between the MID and the 
Gaussian AMID. The latter is defined by restricting the 
optimization to bi-local Gaussian measurements, adopt- 
ing the techniques used for discord in [itI [18| . and we 
provide closed analytical forms for it on special fami- 
lies of two-mode Gaussian states. Section |V] presents 
a thorough comparison between MID, Gaussian AMID 
and Gaussian quantum discord on random two-mode 
Gaussian states, highlighting hierarchical and ordering 
relations between the three measures. Upper and lower 
bounds to the Gaussian AMID at fixed Gaussian discord 
are identified. The states where non-Gaussian measure- 
ments are necessary to reach the minimum in the AMID 
are amply discussed and characterized. In Section IVI[ 
we compare the Gaussian AMID with the Gaussian en- 
tanglement of formation (Gaussian EoF) [2^[2^, show- 
ing on the basis of numerical evidence that the Gaussian 
AMID is always greater or equal than the Gaussian EoF 
for all two-mode Gaussian states, and admits also an 
upper bound at fixed Gaussian EoF (similarly to what 
observed for discord tlTj]). A visual classification of the 
special class of symmetric squeezed thermal states based 
on their degrees of Gaussian EoF, Gaussian discord and 
Gaussian AMID is also provided. Finally, Section IVIII 
summarizes the results we obtained, underlining the main 
implications and delivering hints on possible future ap- 
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plications. 

II. PRELIMINARIES 

A. Gaussian states and measurements 

We consider two modes A and B described by the vec- 
tor — {xa,Pa,xb,Pb)'^ of quadrature operators Xj,pj, 
j — A,B, satisfying the canonical commutation rules 
that can be expressed as [Cj,^fc] = i^jk, where 



il^ J®J, J = 



1 

-1 



(1) 



A two-mode Gaussian state pab can be represented in 
phase space by a Gaussian Wigner function 



W{r) 



1 



TT^ Vdet 7 



(2) 



where r = {xa,Pa,xb,Pb)^ is the radius vector in phase 
space, d = Tt{pabO is the vector of phase-space displace- 
ments and 7 is the covariance matrix (CM) with elements 
7jfe = 2ReTr[pAB(4 - " 4)], j, fc = 1, . . . , 4. Any 

CM 7 has to satisfy the constraints 7 > and 7 -I- iJl > 
(positivity and uncertainty principle) to ensure that it is 
associated to a physical Gaussian state. 

The CM contains a complete information about the 
correlations in a given Gaussian state [l5| . By means 
of local symplectic (unitary on the Hilbert space) oper- 
ations, which leave correlations and entropic quantities 
invariant, the CM of a two-mode Gaussian state can al- 
ways be reduced to a simple standard form 



A C 



/ a ci 

a C2 

ci 6 

V C2 b 



(3) 



where we can assume ci > |c2| > without any loss of 
generality. States with a — b are said to be symmet- 
ric, while states with C2 = ±ci constitute the important 
family of squeezed thermal states. Pure two-mode Gaus- 
sian states are special instances of symmetric squeezed 
thermal states, with zero temperature (det7 — 1), i.e., 
they are locally equivalent to two-mode squeezed vacuum 
states, with standard form covariances 



b — cosh(2r) , ci = — C2 = sinh(2r) 



(4) 



with r being the squeezing parameter. 

Gaussian states can be produced, manipulated and de- 
tected in the laboratory with a high degree of control 
p^ . Some measurements such as photon counting turn 
a Gaussian state into a non-Gaussian one. On the other 
hand, Gaussian measurements play a special role as be- 
ing those that map Gaussian states into Gaussian states. 
These measurements coincide with the standard toolbox 



of linear optics, i.e., can be realized using beam split- 
ters, phase shifters, squeezers, appending auxiliary vac- 
uum states, and performing balanced homodyne detec- 
tion (BHD). Any such a measurement is described by a 
positive operator valued measure (POVM) of the form 



H? 



id,) 
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J 



A.B. 



(5) 



Here the seed element H^ is a normalized density matrix 
of a generally mixed single-mode Gaussian state with 
CM 7j and zero displacements, Dj{dj) = exp(— idjj^) 
stands for the displacement operator, where = 
{xj,pj)^ and dj — {d^'^\ d'^^^ is a vector of certain lin- 
ear combinations of the measurement outcomes of BHDs. 
The POVM ([5]) satisfies the completeness condition 



1 

2^ 



D,{d,)UfD]id,)d'd,^l, 



(6) 



— Ari^^^Ari'-P') following from Schur's lemma 
27| and the normalization condition Tr[n^] = 1. 



where d'^d, — ddf'dd' 



B. Measures of quantum correlations 

Here we briefly review some of the most prominent 
measures of quantum correlations beyond entanglement, 
recently proposed to identify the genuinely nonclassical 
portion of the total correlations in generally mixed states 
of bipartite quantum systems. Achieving a proper un- 
derstanding of the structure and nature of correlations 
in relevant systems is important for gaining insights into 
foundational aspects of quantum theory, and is partic- 
ularly imperative in view of practical applications, as 
quantum correlations yield the key resources to overcome 
classical systems in quantum information protocols. For 
completeness, we also recall the definition of the entan- 
glement of formation. 

Quantum mutual information. — The total amount of 
(classical and quantum) correlations in the state of a 
bipartite quantum system can be reliably quantified in 
terms of the quantum mutual information (28l . [29j 



^qiPAB) = S{pa) +S{pb) - S{pab) , 



(7) 



where S{p) — — Tr(/51n/9) is the von Neumann entropy 
and pa.b are the reduced states of subsystems A and _B, 
respectively. The quantum mutual information Tq (p) of 
a generic Gaussian state p can be easily calculated using 
the formula for the von Neumann entropy of an A'^-mode 
Gaussian state p [30], S{p) = X^ili T?(^i)j where Vi are 
the symplectic eigenvalues ^] of the CM of the state 
and 
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For two-mode Gaussian states PAB^ the two global 
symplectic eigenvalues v±. are defined by 2v'j^ = A ± 

-4det7, with A = det A + det B + 2 det C, see 
Eq. ([3]) [13] ■ The mutual information of a two-mode 
Gaussian state with CM 7 as in Eq. ([3]) is thus 

=i?(VdetA) +;?(x/detB) - 5 - 5 (i^-) , 

(9) 

with S-(x) defined in Eq. ©. 

One-way classical correlations and quantum discord. — 
For any bipartite state whose correlations are purely clas- 
sical, the mutual information can be equivalently ex- 
pressed in two alternative forms 



J^{pab) = S{pa) - inf H^n^^iA\B) , 



{n.} 



J^ipAB) = S{pb) - inf Hrf^MA) , 
{n.} 



(10) 



with 'H^f^^T^{A\B)=J2iPiS{p\^B) being the quan- 
tum conditional entropy associated with the post- 
measurement density matrix p^^^g — TTB[[liPAB]/Pi, 

obtained upon performing the POVM {11^} on system 
B {pi=TT[IliPAB])', the optimization over the POVMs is 
necessary to single out the least disturbing measurement 
to be performed on one subsystem, so that the change 
of entropy on the other subsystem yields a quantifier of 
the correlations between the two parts. 

For arbitrary bipartite quantum states pab , including 
all entangled states and almost all separable states as well 
P], the three quantities in Eqs. ([7]) and pUj) evaluate to 
different results in general, with Iq > J'^, ^7^, and the 

quantities can be interpreted as 'one-way classical cor- 
relation' measures \4^. Such a discrepancy is now recog- 
nized as a signature of nonclassicality of the correlations 
in a given state, and the difference between total correla- 
tion [Eq. ([7])] and one-way classical correlation [Eq. (flUl) ] 
defines what Ollivier and Zurek baptized as the 'quantum 
discord' ji]. 



T^^iPAB) = Iq{pAB) - J^{pab) 

= S{pb) - S{pab) + inf n^f,,M\B) ; 

{n.} ^ 

V^{PAb) - Tq{pAB)-J^{pAB) 



(11) 

(12) 



= S{pa) - S{pab) + inf 'H.f.MA). 

{n,} ^ '-^ 

Quantum discord is an asymmetric measure of quantum 
correlations which has recently received its first oper- 
ational interpretations in terms of the quantum state 
merging protocol [l3 | and has spurred a great body of re- 
search triggered by the investigation of its potential role 
as the key resource for mixed-state quantum computation 
d, 0-11, m, 113] ■ The presence of the optimization over 
local measurements in the definition of discord makes 
its analytical evaluation very hard for general bipartite 
states. No closed formulas are known for TD on arbit rary 
two-qubit mixed states, other than special cases [35|- 



Very recently, a Gaussian version of quantum discord has 
been defined [l3, [3 , where the optimization is restricted 
to Gaussian POVMs of the type Eq. (O, and its closed ex- 
pression has been derived for arbitrary two-mode Gaus- 
sian states [17|. As a consequence of that analysis, it was 
established that the only classically-correlated Gaussian 
states are product states, which are completely uncorre- 
lated. However, in the limit of diverging mean energy 
there can exist two-mode Gaussian states that asyrnp- 
totically approach so-called classical-quantum states [y^, 
where T)'^ — >■ while > [l3l- A symmetrized ver- 
sion of quantum discord — or 'two-way discord' — can 
be defined as 



'D^ipAB) ^Ta&yi{'D^{pAB),'D^{PAB)} , 



(13) 



and in this form becomes vanishing if and only if a state is 
purely classically correlated J^j- Throughout the paper, 
TD will in general denote the Gaussian quantum discord 
[13, Ell, unless explicitly stated. 

Measurement- induced disturbance. — In order to over- 
come the difficulties involved in the evaluation of quan- 
tum discord, Luo introduced the 'measurement-induced 
disturbance' (MID) as an alternative nonclassicality in- 
dicator for bipartite quantum states [l^l- MID is moti- 
vated by the observation that in classical systems, local 
measurements do not induce disturbance. In particular, 
a bipartite state containing no quantum correlations is 
left invariant by the action of any bi-local complete mea- 
surement. On the other hand, even when a state pab 
is a priori nonclassical, any complete bi-local measure- 
ment makes it classical as a result of a decoherence-by- 
measurement process [l3]- MID is thus defined by re- 
stricting the attention to the bi-local complete projective 
measurement £a®£b determined by the eigen-projectors 
£j{k) of the marginal states pj — ^k£j{k) {j = A B), 
where are corresponding eigenvalues, and reads [lO| 



M{pAB)^Ig{pAB)-Iq[£{pAB)], (14) 



where 



£{pab) =Y.P^B{k,l)£A{k) ® £b{1) (15) 

k,l 

is the post-measurement state after local measurements 
£a and £b and pAB{k, I) — TT[pAB£Aik) (8) ^s(0] is the 
probability of obtaining the outcome (kl). The post- 
measurement state is obviously fully classical which im- 
plies that its quantum mutual information ([7]) coincides 
with the classical mutual information of the distribution 
PAB given by I{A : B) = H{pa) + H{pb) - UipAB) M 
with % being the Shannon entropy, where pA and pb are 
reduced probability distributions of the distribution pab ■ 
Hence we can rephrase MID as 



M{pAB)^Iq{pAB)-AA : B) 



(16) 



The MID quantifies the quantumness of correlations in 
terms of the state disturbance after local measurements. 
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but with important differences compared to quantum dis- 
cord: (i) both subsystems are locally probed; (ii) there 
is no optimization over the local measurements, which 
are chosen to be the marginal eigen-projectors for every 
quantum state ^ . Although such a quantity is easily com- 
putable in arbitrary-dimensional systems, and has found 
widespread applications in several investigations [37], a 
number of studies have pointed out that MID is clearly an 
unfaithful and non-refined measure of the nonclassicality 
of correlations in bipartite states [ill, , being nonzero 
and even maximal for states approaching the classical 
limit, and thus severely overestimating quantum corre- 
lations. To date, the MID has not been computed for 
Gaussian states. 

Classical mutual information and ameliorated 
measurement-induced disturbance. — To cure this 
major drawback of MID, one can define an 'amelio- 
rated measurement-induced disturbance' (AMID) by 
incorporating into Eq. (jl6p an optimization (precisely, 
a minimization) over the joint bi-local POVM measure- 
ment Hyi fts on subsystems A and B ll2|. The AMID 
then can be defined as ^, illl, [l^ 

A{pab) = inf {Iq{pAB)-AA: B)} (17) 

= IqipAs) -Ic{pab) , 

where 

2:c(/5ab) = sup I{A : B) (18) 

is the classical mutual information of a quantum state 
PAB where I{A : B) is the classical mutual infor- 
mation of the joint probability distribution pAB{k,l) = 
Tr[/5^BlIyi(fc)®IlB(^)] of outcomes of local measurements 
11^1 and lis on pAB- The AMID captures the quantum- 
ness of bipartite correlations as signaled by the minimal 
state disturbance after optimized local measurements. 
As such, it is a symmetric, strongly faithful nonclassi- 
cality measure jl2| that vanishes if and only if a bipar- 
tite state PAB is genuinely classically correlated @, [ll| , 
and it is operationally interpreted as the quantum com- 
plement to the classical mutual information [Eq. p^ ]. 
while the latter is in turn a bona fide measure of classical 
correlations in general bipartite quantum states [23| . The 
AMID thus englobes the nice properties of discord and 
MID without showing their genetic weaknesses [l^, [3 ■ 
The evaluation and properties of AMID have been re- 
cently investigated for two-qubit systems [11]. ^ 



^ Notice also that this choice of measurements makes MID not 
uniquely defined on bipartite states whose reduced density ma- 
trices have a degenerate spectrum [Till , as it is the case for states 
with maximally mixed marginals. 

2 In Ref. [3 the AMID is defined via an optimization over lo- 
cal projective measurements, rather than more general local 



Hierarchy of nonclassical correlations. — The three en- 
tropic nonclassicality indicators introduced above (dis- 
cord, MID and AMID) satisfy the following hierarchi- 
cal relationship on arbitrary bipartite quantum states 
Pill: 

{V^.V^} <V'^ <A<M. (19) 

Entanglement of Formation. — For pure bipartite 
states \ii)AB, all the measures of nonclassical correla- 
tions introduced above (discord, MID and AMID) reduce 
to the canonical entanglement measure, the 'entropy of 
entanglement', 

E{\'^)ab) = 5(Trs|V'UB(^|) . (20) 

This shows that quantum correlations are faithfully iden- 
tified with just entanglement in the special case of pure 
states of composite quantum systems. 

For bipartite mixed states paBi let us recall that the 
Entanglement of Formation (EoF) Ef{pAB) is a well- 
known entanglement monotone [38| . defined as the con- 
vex roof of the pure-state entropy of entanglement 

EfipAB)^^ niin Vp.i?(|^)jiB) , (21) 

where the minimum is taken over all the pure-state real- 
izations of PAB , 

PAB = J^P'^'^^abWab ■ 

i 

For general mixed bipartite Gaussian states pab, one 
can introduce a Gaussian version of the EoF (Gaussian 
EoF) E'^ , which is by construction an upper bound to 
Ef, defined as the convex roof of the entropy of entan- 
glement restricted to decompositions of pab into pure 
Gaussian states only [2^. The Gaussian EOF can be 
evaluated via a minimization over CMs: 

Efh)^ inf E{Y), (22) 

7 <7 : dct(7') = l 

where the infimum runs over all pure bipartite Gaussian 
states with CM 7' smaller than 7. Compact formulas 
for E'^ exist for all symmetric two-mode states (where 
i?^ ^ Ej: as the Gaussian decomposition is globally op- 
timal) (25} ■ and in the nonsymmetric case for squeezed 
thermal states and so-called states of partial minimum 
uncertainty [3^; for all the other two-mode states its 
value can still be found analytically '26^ . 

In the following Sections, we will approach the evalua- 
tion and characterization of MID and AMID for arbitrary 
two-mode Gaussian states and their interplay with en- 
tanglement, complementing the studies of Refs. [13, 13 
on discord. 



POVMs. Both versions of AMID have also been studied in 
Ref. , albeit without naming the considered measures explic- 
itly 
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III. MEASUREMENT-INDUCED 
DISTURBANCE OF TWO-MODE GAUSSIAN 
STATES 

In this Section we will address the calculation of the 
MID M [3, Eq. (HH), for two-mode Gaussian states. 
Here and in the following we will always consider, with- 
out any loss of generality, states pab whose CM 7 is in 
standard form, Eq. ([3|). The expressions we derive can 
be straightforwardly recast in terms of a set of four local 
symplectic invariants for general two-mode states, that 
uniquely define the standard form covariances (32l . [40| . 

The quantum mutual information of a two-mode 
Gaussian state can be computed via the formula 
Here we will be concerned with the calculation of the 
mutual information after local projections onto the eigen- 
states of the marginal density matrices. 

For a generic two-mode mixed Gaussian state pab in 
standard form, the reduced states are just thermal states 
/5th, A and pth,_B with mean number of thermal photons 
("a) = (a — l)/2 and (ns) = {h — l)/2, respectively. 
The local measurements £a <8> £3 entering the expression 
of MID, Eq. (|14l) . are then non-Gaussian measurements, 
given by projections onto Fock states (joint photon count- 
ing), 



and the post-measurement state reads as 



£{pab) = ^ p{m,n)\ 



mjAiml (g) \n)B{n\ 



(23) 



(24) 



where p(m, n) =A{'m'\B{n\pAB\'m') aI^) b is the joint prob- 
ability distribution of finding m photons in mode A and 
n photons in mode B. The determination of the quantum 
mutual information l!g[£{pAB)] of the post-measurement 
state (IM)) requires the calculation of the local von Neu- 
mann entropies S{p^), j — A,B of the reduced states 

Pa,b 



= Tr 



b,a 



£{pab) 



as well as the global entropy 

S[£{pab)]- Since the reduced states are just equal to the 
local thermal states, i.e. /5| = Pthj", j — A, B, we obtain 
immediately the local entropies equal to S{pa) — 
and 5(/5|) = ^J(6), where ^{x) is defined by Eq. dH). The 
global entropy can be computed from the joint photon- 
number probability distribution p{m,n). The latter can 
be derived using the generating function for the distribu- 
tion in the spirit of Ref. 14H. 

We start by noting that any two-mode quantum state 
PAB can be represented by the complex normal quantum 
characteristic function defined as 



C(/3i,/32) = TiipABef'^'+f^-^'e-^'^^-f''-'''] 



(25) 



where a (d^) and b (b^) are annihilation (creation) oper- 
ators of modes A and B, and /3i, are complex param- 
eters of the characteristic function. For a Gaussian state 



in the standard form ^ with zero means (d) = (&) = 
the characteristic function attains the form 



C(/3i,/?2) = exp[-(Bi|/3i|2 -1-^21/321') 
+Dpip; + c.c.)] , 



(26) 



where Bi = (AdtAa) = (a - l)/2, B2 = (AS^AS) = 
(6- l)/2, D = (AaA6) = (ci -C2)/4, D = -{Aa^Ab) = 
— (ci + C2)/4. Here c.c. stands for complex conjugate 
terms and AA = A — (A) . The generating function for 
the distribution p(m, n) then reads [4^ 



G(Ai,A2) 



1 



7r2AiA2 



exp 



M 

Ai 



1/32 P 



C(/3i,/32)d2/3id2/32 



(27) 



where Ai and A2 are real parameters. By inserting the 
characteristic function (|26p into the integral (|27|) and per- 
forming the integration, we arrive at the generating func- 
tion G(Ai,A2) = Fi (Ai,A2)^2 (Ai,A2), with 



i^,(Ai,A2) = 



^l + BiAi+B2A2 + ifjAiA2' 



(28) 



where = \[a - 1)(6 - 1) - c2]/4, j = 1, 2. The sought 
photon-number distribution is obtained by differentiating 
the generating function ((Tf)) as 



p(m, n) 



(-l)"+"9™+"G(Ai,A2) 



7 In I 



(29) 



Ai=A2=l 



which gives us explicitly 

m\n\ ^-^ ^-^ \v\l \V2} 

1/1=01/2=0 \ \ ^/ 

Y.Q^^\m-v^,n^V2), (30) 



where 

g(^)(a,/3) = 



(i3i+g,-)"(i?2+j^,)^ 

4a+/3(l + B1+B2+ Kj)°'+I^+^ 

mi^,,) [2(a + /3-/)] 

2. ia + l3-l)l 



1=0 

1 + B,+B2+K, 
\B,+K,){B2 + K,) 



(31) 



The global posterior von Neumann entropy is then given 
by the formula 



S[£{pab)] = - ^ p{m,n)\np{m,'i 



(32) 



m,n=0 



and can be evaluated numerically for arbitrary two-mode 
Gaussian states. 
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We note that in Refs. [43| it is shown that the 
joint photon-number distribution for multimode Gaus- 
sian states can be written in terms of multivariable Her- 
mite polynomials. Recursive formulas can be derived for 
the calculation of the multivariable Hermite Polynomials, 
and this fact can be exploited to speed up the numerical 
evaluation of the MID. 

The final expression for the MID of two-mode Gaus- 
sian states pAB with a standard form CM as in Eq. ([3]) 
is then 

M{pab) = S[£{pab)]-S{pab) (33) 
= - ^ p{m,n)\i\p{rri,n) - - ^{v-) . 

m.n—O 

In the subcase of two-mode squeezed states, with ci = 
±C2 = c, we have Ki = K2 ^ K = [{a - - 1) - c^]/4: 
and the distribution ([50)) takes the simplified form 



p{m,n) = 



{Bi+Kr{B2 + KT 
m\n\{l + B1+B2 + 
min(m,n) 

3=0 



1 + Si 



3 

B2 



j\{m + n- j)\ 
+ K 



iBi+K){B2+K)\ ■ ^^^^ 

For the very special instance of pure two-mode Gaus- 
sian states with the standard form given by two-mode 
squeezed vacuum states 



(35) 



71=0 



with standard form covariances given by Eq. Q, where 
q — tanhr, the MID can be evaluated in closed form. 
The state is in Schmidt decomposition with Schmidt 
coefficients A„ ~ ^1 — g^g". Being a pure state, 
S{\ip{r))AB) — 0. The post- measurement state then 
reads [ioj 



£imr))AB) = il-cf)J2q' 

n=0 



\n)A{n\(E)\n)B{n\, (36) 



and its von Neumann entropy is S[£{\'ip{r)) ab)] — 
— ^„A^jlnA^j — 5^[cosh(2r)], which precisely coincides 
with the entropy of entanglement of the pure two-mode 
Gaussian state, Eq. I^U^. Therefore, 

MU{r))AB] = Emr))AB]^ cosher) ln[cosh\r)] 
— sinh^(r) ln[sinh^(r)] , (37) 

as expected from the definition of MID. Let us remark 
that this value is attained by local non-Gaussian mea- 
surements (joint photon counting). An interesting ques- 
tion is whether there exist local Gaussian POVMs that 
can result in a measurement-induced disturbance equal to 
the one in Eq. p7p for pure two-mode Gaussian states. 
Surprisingly enough, we will prove in the next Section 
that the answer is negative. 



IV. GAUSSIAN AMELIORATED 
MEASUREMENT-INDUCED DISTURBANCE OF 
TWO-MODE GAUSSIAN STATES 

Given the crucial role of Gaussian states in quantum 
information processing, it is important to carry out a 
thorough analysis and comprehensive characterization of 
their quantum and classical correlations. A primal step 
to undertake is to approach the evaluation of faithful 
measures, such as the classical mutual information Tc and 
correspondingly the AMID A [Eqs. (|17ll8p ]. for arbitrary 
two-mode Gaussian states pab- However, the problem 
appears formidable if any possible non-Gaussian POVM 
is allowed. Therefore, similarly to what has been done 
for discord [13, [3 j one can define a Gaussian version of 
the AMID, A'^ , where the optimization is constrained to 
local Gaussian POVMs of the form ([5]), as 



A'^ipAB) ^MpAB)~:i:ciPAB), 



where 



IciPAB) 



sup I{A : B) 



(38) 



(39) 



is the Gaussian classical mutual information of the quan- 
tum state Pab- The true AMID, optimized over general 
local measurements, would be then bounded from above 
as 



A{pab) < niin {A'-{pAB),MipAB)}- 



(40) 



A relevant and natural question, related to the one in 
the closing of the previous Section, is whether Gaussian 
measurements are always optimal for the evaluation of 
the AMID, i.e., whether A = A^ for all bipartite Gaus- 
sian states. In the case of discord, there are subclasses 
of two-mode Gaussian states where one can prove that 
the Gaussian discord achieves the global minimum in 
Eq. (jlip even when including potentially non-Gaussian 
measurements on a local subsystem 17[: the properties 
of the Gaussian discord allow us to conjecture that this 
might be true for arbitrary two-mode Gaussian states, 
although no rigorous investigation is available to date 
to support this claim ^ . Remarkably, we will find in- 
stead that when both parties are locally probed, there is 
a finite- volume set of two-mode Gaussian states (notably 
including pure states) for which no n- Gaussian measure- 
ments are necessarily optimal for the evaluation of the 
AMID, and the Gaussian AMID only provides a strict 
upper bound to it. This proves that Gaussian joint mea- 
surements may not be the least disturbing ones on general 
two-mode Gaussian states. 



^ A preliminary numerical study reveals that, for general two- 
mode mixed Gaussian states, even allowing for a non-Gaussian 
measurement such as photon counting, the corresponding discord 
is found to be never smaller than the one associated to the opti- 
mal Gaussian measurement [A. Datta, private communication]. 
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We hereby develop the framework for the determina- 
tion of the Gaussian AMID on general two-mode Gaus- 
sian states, and provide closed formulas for it in some 
special cases. The nontrivial part in the determination 
of the quantity (pS]) is the calculation of the classical mu- 
tual information I^{pab) requiring maximization of the 
Shannon mutual information X{A : B) over local Gaus- 
sian POVMs i\A and Hs of the form 

We begin by observing that in the considered optimiza- 
tion task we can restrict ourselves to covariant Gaus- 
sian POVMs ([5]) projecting onto pure states (rank-one 
POVMs), similarly to discrete- variable scenarios 
The proof of this statement is provided in Appendix 
Rl We can thus focus on local measurements of the 
form ([SJ, where the seed state Ilj is a pure single- 
mode Gaussian state with CM 7^. Let us recall that 
any pure-state one-mode CM 7^ can be expressed as 
7^- = U{ej)V{rj)U'^{ej), where 





-2ri 



(41) 

where 9j S (0, tt) and rj > 0. In this picture, homodyne 
detection on mode j is recovered in the limit of an in- 
finitely squeezed pure state Ilj, i.e., rj — >■ 00. On the 
other hand, heterodyne detection on mode j corresponds 
to rj = 0. 

Now we want to maximize the Shannon mutual infor- 
mation X{A : B) of the distribution (IA3p . 

P{d) = TliflAidA) <E} UB{dB)pAB], 

over all single-mode pure-state CMs 7a, s- Expressing 
the two-mode state CM 7 in block form as in Eq. (jSj, and 
using the formula for the Shannon entropy of a Gaussian 
distribution P oi N variables with classical correlation 
matrix S, HiP) = ln[(27re)f Vdet S] [s^, the sought 
mutual information can be obtained in the form 44 1 : 



I{A:B) = ^ln 



det {jA + A) det {-/b + B) 
det {jA ® 7b + 7) 



(42) 



Since the determinant is invariant with respect to sym- 
plectic transformations we can assume the CM 7 to be in 
standard form, Eq. ([3]); moreover, given Eq. ([8]), Eq. (|411) . 
the invariance of the determinant under orthogonal trans- 
formations, and the monotonicity of the logarithmic func- 
tion, the object to be maximized reads 

^ ^ ^ det A' det 

f{rA,rB,dA,eB) = J— , 43 

det 7' 

where the CM 7' has the 2x2 blocks of the form 
A' ^al + V{rA), B' = bl + V{rB), 
C ^ U^ieA)diag{ci,C2)UieB) . 



(44) 



We recall that the Gaussian AMID takes then the form 
with 



1 



(/Sab) = 77 In 



sup firA,rB,t 



(45) 



The determinant in the denominator of Eq. (|43)) can 
be expressed in terms of the local symplectic invari- 
ants /i = det A', I2 = detB', I3 = det C" and I4 = 
Tv{A'JC'JB'JC'^J) as det 7' = Iih + /| - /4 |i01. Ev- 
idently, the function p3)) depends on the phases Oa,b 
only through invariant li so we can optimize / over the 
phases by optimizing I 4^. Since I 4 is always nonnegative, 
the maximum of the function / will be obtained if the 
phases will maximize I^. 

Three different cases must be distinguished depending 
on the values of the eigenvalues ci,2 of the matrix C. 



1. If ci — C2 — 0, then 13.4 = 0, which implies 
f — \ and therefore I^{pab) — 0. Since also 
^qiPAB) = 0, we get finally A'^Ipab) = 0. Ob- 
viously, product states have neither quantum nor 
classical correlations. Actually, it is known that any 
non-product two-mode Gaussian state has nonzero 
quantum correlations, so we must expect be > 
as soon as one of the covariances in the block C 
is nonzero. 

2. If ci > and C2 — 0, then for rA,B > the optimal 
phases are 0A,B,opt = t^/^ and for these phases the 
invariant I4 takes the value I'^ = cl{a + e'^'^^){h + 
e^^"). Let us define 



1 - h{rA,rB) 



= 9{rA,rB) 
= f{rA,rB,( 



opt, 6*5, opt)- (46) 



We have h = I'^/{Iil2) and the function h (and 
consequently g) is obviously maximized in the lim- 
its rA,B — > 00 (doubly homodyne detection) when 
we get 



ah 



Sho 



ab 



(47) 



Since ghom is also larger than the maximal values of 
the function g on the boundaries = or = 0, 
we arrive at the conclusion that (j47p is the optimal 
value of Eq. (|43p for states with C2 = 



3. In the general case Ci,C2 7^ 0, the invariant 
I4 can be expressed as I4 — c\c^Tv{XB'), 
where X is a real symmetric positive- 
semidefinite matrix of the form: X — 
U]jdia.g{c^^,c^^)UAA'U'^d\a.g{c^^,c^^)UB, with 
Uj = U{9j) defined in Eq. (j4T|) . and the diagonal 
matrices A' and B' defined in Eq. (H^ . Expressing 
the matrix X through eigenvalue decomposition 
X = W{(j))dia.g[xii9A),X2ieA)]W^{(j)), where W 



^ Recall that in our convention ci > |c2|, therefore the optimal 
Gaussian measurement consists in the local homodyne detection 
of the quadrature with the highest intermodal correlation. 
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is an orthogonal matrix diagonalizing X and 
Xi{Oa) > X2{0a) are the eigenvalues of X depend- 
ing on the angle 9a, we can maximize /4 over 
the phase (j). Further, if we substitute into the 
obtained formula the explicit forms of eigenvalues 
xi.2{0a), we can perform also the maximization 
over the angle 9 a, which finally yields the invari- 
ant /4 maximized over the phases 9a,b of local 
measurements of the form: 



4 = {[a + cosh(2rA)] [6 + cosh(2rB)] 
-I- sinh {2rA) sinh {^tb)} (c? + c^) 
+ {[a + cosh {2rA)] sinh (2rB) 
+ [b + cosh (2rB)] sinh {2rA)} {c\ 



(48) 



The corresponding value of / is denoted by 
g{rA, tb)- For a generic two-mode mixed Gaussian 
state it is again convenient to express the function 
gifA, Tb) in terms of h{rA, tb) as in Eq. being 



where 



h{rA,rB) = {Ii-ll)l{hh), 



h = (a + e2'^-^)(a + e-2-^), 



h = {b + e 

h = ClC2 



-2rB 



(49) 



(50) 
(51) 
(52) 



Since h > as follows from the inequality 2{A : 
B) > 0, we have to maximize h. Introducing new 
variables A = e^*"^ and fi = e^''^, the extremal 
points of h can be found by solving the stationarity 
conditions dh/dX = and dh/dfi ~ 0, respectively, 
leading to a set of coupled polynomial equations of 
the form: 

cfia + A)V^ + [clb{a + - c^6(aA + 1)^ 
+clcla{X^ - - cl{aX + 1)^ = 0, 
cl{b + fifX^ + [cla{b + fif - cjaibfi + i f 



+clclb{i? - 1)]A - clib^i + if 



0. 



(53) 



Upon solving the first equation as a quadratic equa- 
tion with respect to /x and substituting the obtained 
roots into the second equation, one arrives after 
some algebra at a single 12"'-order polynomial in 
the variable A that we do not write here explicitly 
due to its complexity. By taking its real roots cal- 
culated numerically together with stationary points 
on the boundary and picking the one for which h is 
maximized, we can finally get the optimal squeezing 
parameters ta^b of the seed elements Ha.b of opti- 
mal local POVMs [Eq. ([5])] maximizing the classical 
mutual information and thus attaining the Gaus- 
sian AMID, Eq. (1551) . of a generic two-mode Gaus- 
sian state. 



Analytical progress in the calculation of can be 
achieved for special classes of two-mode states, detailed 
in the following. 

4. Symmetric states. For symmetric states with a — 
b the maximization of can be in principle per- 
formed analytically. As the optimal solution is 
clearly symmetric with ta — tb = r, we have 
only one stationarity condition dh/dX = 0, where 
A = e^'' is the parameter to optimize over. After 
some algebra the stationarity condition boils down 
to the following fourth-order polynomial equation 



a^X -I- as A -t- 02 A +aiX + aQ — 0, 



(54) 



where 

ao = 

04 — 
03 



2 { 1 ^ 

—C2, ai ~ —a (Cj^Cr 

2 o 2 / 2 

Ci, 0,2 — Jfl (Cx — t 

22,02 2 2\ 
' 3ci - a C2J 



3c2 - a^c?) 



= a CiC 



l'-2 



(55) 



The equation (IM|) can be solved analytically using 
the Cardan formulae but the obtained solutions are 
rather cumbersome and therefore we do not give 
them here explicitly. Calculating the values of the 
function h{X) in the admissible real solutions of 
Eq. ([5^ and also in the stationary points on the 
boundary, the point in which the function is maxi- 
mal gives us the sought optimal squeezing. 

5. Squeezed thermal states. For generally nonsym- 
metric squeezed thermal states with ci = ±C2 = c, 
the optimal Gaussian POVMs can be derived in a 
simple closed form by performing the maximization 
of h in Eq. p9)) . We find the following results. If 
the state parameters a, b and c satisfy the inequality 
(a -I- & -1-1)^ > ab{ab — c'^), then the optimality is ob- 
tained by homodyne detection (r^^s 00) on both 
modes giving Eq. (|47)) : conversely, if {a + b + lf < 
abiab— c^), then the optimality is obtained by het- 
erodyne detection [rA = = 0, projection onto 
coherent states) on both modes giving 



5het 



(a + l)(& + l) 



(a + l)(fe+l) -c2 



(56) 



Summarizing, the Gaussian AMID of two-mode 
squeezed thermal states is given by Eq. (|45p with 



.ghom [Eq. dUl) 

ffhct [Eq. m] 



[a + b+lf > abiab-c^); 
otherwise. 



(57) 



Pure states.. For pure two-mode squeezed vacuum 
states with a ~ b ~ cosh(2r'),ci — — C2 = c — 
sinh(2r), the first case in Eq. (|57p is always satis- 
fied, and doubly homodyne measurements are opti- 
mal for the calculation of the Gaussian AMID. Al- 
though this result is quite intuitive (one could guess 
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that in the pure-state case the optimal local Gaus- 
sian measurements possessing maximum Shannon 
mutual information for distribution of their out- 
comes would be homodyne detections), the corre- 
sponding value of A'^ is strictly bigger than the en- 
tropy of entanglement, Eq. (|20l). which corresponds 
to the true AMID A globally optimized over joint, 
possibly non-Gaussian local POVMs as in the def- 
inition (ITT)) : and we know from Eq. (|37|) that the 
latter is indeed attained by local photon counting: 
A^mr))AB) > A{mr))AB) = X(IVW)ab) = 
E{\i;{r))AB)- Namely, 

^'^(I^W)ab)) = 2Mmr))AB) - ln[cosh(2r)] , (58) 

which is strictly bigger than the expression in 
Eq. dSll). 



The latter finding entails a novel and interesting show- 
case in which non-Gaussian local measurements on Gaus- 
sian states can lead to the extraction of larger corre- 
lations than any pair of local Gaussian measurements 
if quantified by mutual information. In the case of 
pure two-mode squeezed states this behavior cannot be 
simply attributed to the fact that, while the state en- 
codes perfect correlations between photon numbers (as 
expressed by Pn,m = for n ^ m), there are imper- 
fect correlations between quadrature operators (as ex- 
pressed by a nonzero Einstein-Podolsky-Rosen variance 
{{xA —xbY) = e~^^). This is because, while such a vari- 
ance decays with increasing squeezing r, the gap between 
the Gaussian AMID A'~^ (obtained upon local homodyne 
detections) and the true AMID A = A4 (obtained upon 
local Fock projections), for pure states, is a monotoni- 
cally increasing function of r converging to 1 — In 2 w 0.3 
in the limit r — ^ cxd. 

A full numerical comparison between Gaussian AMID, 
MID and Gaussian discord for two-mode Gaussian states 
will be provided in the next Section. 



V. COMPARISON BETWEEN 
NONCLASSICALITY MEASURES FOR 
GAUSSIAN STATES 

Here we deploy a comprehensive comparative analy- 
sis of MID M [Eq. (HH), Eq. Gaussian AMID A^ 
[Eq. ([HS]). Eq. (US])], and two-way Gaussian quantum 
discord [Eq. Ref. [13], as tools to quantify 
the quantumness of correlations in arbitrary two-mode 
Gaussian states via entropic descriptions of the state dis- 
turbance following suitable local measurements on one 
or both local parties. The results of the previous Sec- 
tions show that in Ineq. (j40p . either quantity on the right 
hand side can be the smallest on particular instances of 
two-mode Gaussian states, suggesting that the subset of 
two-mode Gaussian states whose true AMID A, Eq. ([T7| , 
is necessarily optimized by non-Gaussian measurements, 
might have a finite volume in the space of general two- 
mode Gaussian states. 



To confirm this interesting feature and to investigate 
the tightness of the hierarchy established in Eq. p^ . 
we have generated a large number of random two-mode 
Gaussian states (up to 10^), and for each of them we 
have evaluated the three symmetric nonclassicality indi- 
cators 2?^, A^ , and A^, following the prescriptions of the 
preceding Sections. The resulting analysis is illustrated 
in Fig. [1] Panel (a) shows that, while the MID can be 
arbitrarily larger than the Gaussian AMID in principle, 
there is nonetheless a finite region in the {A'^ , A4) dia- 
gram that allocates two-mode Gaussian states for which 
even non-optimized non-Gaussian measurements (specif- 
ically, photon counting) result in a larger classical mu- 
tual information, hence minimize the quantum correla- 
tions in the definition of the AMID, compared to the 
optimal Gaussian POVMs. In this study we did not 
consider other non-Gaussian measurements apart from 
photon counting (that enters in the definition of MID 
because Gaussian mixed states are essentially thermal 
states) : Therefore, we can expect that the region in which 
general non-Gaussian measurements are optimal for the 
AMID can be in principle much larger than the one high- 
lighted by the present study (that is located below the 
blue online segment of equation A'-^ = in Fig. 1(a) I. 



Still, our finding is perhaps one of the most striking in- 
stances of an operational quantum informational mea- 
sure for Gaussian states that can gain a significant opti- 
mization by the use of suitable non-Gaussian operations. 
Non-Gaussian operations can sometimes reveal quantum- 
ness more accurately, thus unleashing more precisely the 
available nonclassical resources, than the best Gaussian 
measurements, on certain two-mode Gaussian states, in- 
cluding quite remarkably all two-mode pure Gaussian 
states. Notice however that the rigidity in choosing the 
non-Gaussian measurements for the evaluation of MID, 
excluding any optimization procedure, results in most of 
the cases into a very loose overestimation of the quantum 
correlations, as testified by the unbounded region above 
the straight blue line, filled by states where certainly joint 
photon counting is not optimal for the AMID. Interest- 
ingly, pure sta tes em body the lower bound (dashed black 
curve) in Fig. l(a)[ they are therefore the states where 
the Gaussian AMID realizes the most dramatic overes- 
timation of the true AMID, that nonetheless can never 
exceed sa 0.3 as computed in the previous Section. A 



family of states sitting on the blue line in Fig. 1(a) will 
be characterized shortly. 

Before that, let us comment on the other panels of 
Fig. [TJ Panel (b) shows as expected, and in full analogy 
with the case of two qubits 12], that in general the un- 
optimized MID based on photon counting is a very loose 
upper bound to quantum discord for two-mode Gaus- 
sian states (reducing to it on pure states, depicted as 
dashed black again), unbounded from above and relent- 
lessly approaching arbitrarily large values even for states 
with nearly vanishing quantum correlations as quanti- 
fied by the (Gaussian) discord. This should discourage 
the usage of MID in general as it almost always provides 
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FIG. 1: (Color online) Comparison between (a) MID versus Gaussian AMID, (b) MID versus two-way Gaussian discord, and 
(c) Gaussian AMID versus two-way Gaussian discord, for IC randomly generated mixed two-mode Gaussian states. Pure 
two-mode squeezed states are accommodated on the dashed black curve in all the plots. See text for details of the other 
boundaries. All the quantities plotted are dimensionless. 




FIG. 2: (Color online) Comparison of different measures of quantum correlations for two-mode symmetric squeezed thermal 
Gaussian states (b — a,ci — —C2 = c). Panels (a)-(b): Two-way Gaussian quantum discord "V^ = max{0^,I?~*} (dotted 
black line), Gaussian AMID A'^ associated to optimal bi- local Gaussian POVMs (dashed blue line), and unoptimized MID A4 
associated to joint photon counting (solid red line), plotted versus the normalized state covariance parameter c/V — 1, for 
(a) a — 1.05 and (b) a — 2. The AMID A optimized over all possible (Gaussian and non-Gaussian) measurements is certainly 
A < mm{M,A'^}. Both MID and Gaussian AMID majorize the Gaussian discord, but for c bigger than a certain threshold value 
c*(a) (ticked by a vertical gray line in the plots) one has M < A'^ , meaning that non-Gaussian measurements become necessarily 
optimal for the AMID. Panel (c) depicts the threshold curve c*(a) (solid black line), defined by the condition j\4 = A'^ , in the 
normalized parameter space {{a — l)/a, c/Va^ — 1}. The shaded (orange) region above the threshold line allocates instances 
of the considered family of states, lying in the neighborhood of pure two-mode squeezed states, where certainly Gaussian 
POVMs are not globally optimal for the AMID, since photon counting results in a lower figure of merit. Below the threshold, 
either Gaussian measurements are optimal or there may exist some more general non-Gaussian measurement that achieves the 
absolute minimum in A: our analysis cannot rule out this possibility. All the quantities plotted are dimensionless. 



overestimations, rather than reliable quantifications, of 
nonclassicality of bipartite correlations. 

The last panel shows also a somehow analogous situ- 
ation to the two-qubit case [l^l: the Gaussian AMID is 
intimately related to discord, and admits upper and lower 
bounds at a given value of the two-way Gaussian discord. 
The lower (blue online) boundary in panel (c) accommo- 
dates states for which the two quantifiers give identical 
prescriptions for measuring quantum correlations. These 
are states with CM in standard form, Eq. ([3]), given by 



a = cosh(2s), b = cosh^(r) cosh(2s) -I- sinh^(r) 
Ci = — C2 = coshr sinh(2s) , 



(59) 



in the limit r oo. They are characterized by A'^ — 
= 2sinh^ sln(coths). Pure states fill once more the 
dashed black curve, for which = £, Eq. (|37p . but 
A*-^ is strictly big ger, E q. ([55]) . The upper (red online) 
boundary in Fig. 1(c) can be spanned for instance by 
symmetric squeezed thermal states (b — a, ci — ~C2 — 
c), with a ^ 1 and c € [0, Vo? — 1). Upper and lower 
boundaries ideally conjoin asymptotically for diverging 
discord and Gaussian AMID. 

We can now analyze in detail the competition between 
the MID associated to photon counting (typically very 
loose, but optimal on pure states) and the Gaussian 
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FIG. 3: (Color online) Plot of MID (solid red line), Gaussian 
AMID (dashed blue line) and Gaussian two-way quantum dis- 
cord (dotted black line) as a function of the parameter a for 
the two-mode Gaussian states of Eq. (|60[) . These curves give 
origin to the dotted magenta lines in Fig. [IIa),(b). All the 
quantities plotted are dimensionless. 
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FIG. 4: (Color online) Plot of Gaussian AMID A° versus 
Gaussian EoF Ef for 10^ random two-mode Gaussian states. 
The dashed line of equation A'^ = Ef stands as a lower bound 
for the physically admitted region. Refer to the main text for 
details of the other curves. All the quantities plotted are 
dimensionless. 
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AMID (very accurate for mixed and strongly correlated 
Gaussian states) to maximize the classical mutual infor- 
mation, hence minimizing the AMID, Eq. (|40l) . on two- 
mode Gaussian states. We believe it be relevant to focus 
on the class of two-mode symmetric squeezed thermal 
Gaussian states {b = a,ci = — C2 = c), for which the 
involved measures can be simply evaluated ^. We plot 
in Fig. [2] [(a), (b)] a comparison of the three quantumness 
measures studied in this work as a function of the rescaled 
state parameters a and c. We see that there is a certain 
threshold value c*(a) beyond which the Gaussian POVMs 
are no longer optimal for the AMID, and non-Gaussian 
measurements such as photon counting (via MID) pro- 
vide a more accurate result, culminating in the extreme 
case of pure states where those specific measurements are 
globally optimal. Panel (c) depicts the threshold in the 
parameter space, highlighting the region where our anal- 
ysis conclusively reveals the necessity of non-Gaussian 
measurements for the global optimization of AMID and 
classical mutual information of the considered class of 
two-mode Gaussian states. As previously remarked, this 
region can be in principle (and is likely to be so) much 
larger. Yet, it certainly occupies a finite volume in the 
space of general two-mode Gaussian states. Notice that 
for all the Gaussian states in such a region, non-Gaussian 
measurements allow one to extract stronger correlated 
measurement records compared to any bi-local Gaussian 
measurement, as the classical mutual information is max- 
imized by non-Gaussian detections. The states attaining 
the threshold identified in this analysis, are an instance 



of states filling up the blue line in Fig. 1(a) 



Finally, we exhibit an example family of two-mode 



Gaussian states where, on the opposite end, the MID 
based on non-Gaussian detections is a highly inaccurate 
measure of quantum correlations. These states sit on the 
quasi- vertical dashed (magenta online) curves in Fig.[TJa) 
and (b). They are symmetric states with 



a, C2 = 0, and ci = (a^ — 1 — In a)/a. 



(60) 



^ For squeezed thermal states, the Gaussian discord is optimized 
by a local heterodyne detection [TtI . [Tsll . 



As apparent from Fig. [3J their Gaussian discord and 
Gaussian AMID stay limited (smaller than w 0.06) and 
rigorously vanish in the asymptotic limit a — oo. On 
the other hand, their MID arising from Fock projections 
increases arbitrarily and diverges for a — > cx), embodying 
an extreme overestimation of some vanishing quantum 
correlations. Clearly there will be many more families of 
Gaussian states where such a behavior will arise. 



VI. NONCLASSICALITY VERSUS 
ENTANGLEMENT 

Here we present a numerical comparison between non- 
classicality of correlations, measured by means of the 
Gaussian AMID [Eq. (|38|) ]. and entanglement, quantified 
by the Gaussian EoF [Eq. d^^ ] ,26j, for generally mixed 
two-mode Gaussian states. A similar analysis was per- 
formed in Ref. fv?], with (Gaussian) discord used as a 
nonclassicality indicator. 

Fig. m shows the distribution of Gaussian AMID ver- 
sus Gaussian EoF for a sample of 10^ randomly generated 
two-mode Gaussian states. In analogy with the case of 
Gaussian discord vs Gaussian EoF, it is possible to iden- 
tify upper and lower bounds on the Gaussian AMID A'^ 
at fixed entanglement E'j^. Interestingly, our numerical 
exploration shows that for all two-mode Gaussian states 
PAB, it is 



(61) 
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This provides a novel hierarchical relationship between 
different types of nonclassical resources, entanglement 
E'^ , and more general measurement-induced quantum 
correlations A'^: On the basis of the hereby employed 
measures (both symmetric by construction and restricted 
to a fully Gaussian scenario), the latters appear to always 
encompass and exceed entanglement itself for two-mode 
generally mixed Gaussian states. A similar relationship 
does not hold for discord, which can be smaller as well as 
larger than entanglement of formation, even in a Gaus- 
sian scenario 17, 48]. 

We can provide two families of two-mode Gaussian 
states for which Eq. (|6ip becomes asymptotically tight. 
One such class is provided, e.g., by symmetric squeezed 
thermal states, whose standard form CM is as in Eq. ([3]) 
with 



Cl 



-C2 



where 



j>>0, a > max{i>, (1 + i>2)/(2;>)} 



(62) 



(63) 



The Gaussian EoF of these states (equal to the true EoF 
minimized over all possible pure-state decompositions, 
by virtue of the symmetry of the states [13 ) is a simple 
monotonically decreasing function of the positive param- 
eter i>, 



(1 



(1 



(1- 



4z> 



(64) 

if < 1, and Ef{i>) = otherwise. The Gaussian AMID 
can be computed analytically according to the prescrip- 
tion of Sec. IIV| and for any fixed (i.e., fixed Gaussian 
EoF) one can find the optimal value of the parameter 
a, in the range defined by Eq. (j63l) . that minimizes 
The resulting A'~^ as a function of E'^ is plotted in Fig. |4] 
as a green curve: For this family of states, the Gaussian 
AMID approaches the Gaussian EoF as the latter tends 
to zero, tending to saturate Ineq. (|6ip asymptotically in 
the regime of infinitesimal correlations. 

Furthermore, let us consider another class of symmet- 
ric two-mode states, whose standard form CM is as in 
Eq. §1 with 

b = a,ci^ a-{l+i)'^)/{2a), ca = a~{2a)/{l+i)'^) , (65) 



where the parameter range is the same as in Eq. (|63p. 
The Gaussian EoF of these states is still given by 
Eq. (|64p . while in the limit a — > oo one can show that 
their Gaussian AMID is attained by homodyne detec- 
tions on both modes, yielding A'~^{iy) = 1 — ln(4z>) + 
ln(l -t- 1>^). The corresponding Gaussian AMID vs Gaus- 
sian EoF curve for these states is depicted in Fig. |4]as a 
blue curve: In this case, the Gaussian AMID approaches 
the Gaussian EoF as the latter tends to infinity, also sat- 
urating Ineq. (|6T|) asymptotically. 

It has to be underlined that the two presented Gaus- 
sian families are just examples to show that the bound 



in Eq. ((6T|) can be asymptotically tight, and we remark 
that the combination of the two presented curves does not 
provide a strict lower bound to Gaussian AMID against 
Gaussian EoF for all two-mode Gaussian states, as it is 
evident from the presence of some random points below 
the intersection of the two (green and blue) curves - how- 
ever above the dashed line corresponding to the bound 
of Eq. dnH) [see Fig.|l]. 

On the other hand, a tight upper bound on the Gaus- 
sian AMID at fixed Gaussian EOF can be identified. We 
found it numerically to be constituted by the maximum 
of two branches, 

1 + 21n(l + £>) - ln(4£>), 
4/c - 1 < < 1; 
A^ipAB) < { (66) 
ln(l + t>) - ln(i>), 
< i>< 4/e- 1; 

where E'^{pab) is given by Eq. (pi)) . 

The first expression in Eq. (|66|) corresponds to the 
Gaussian AMID of the states of Eq. ([M]) with s ^ oo 
and r = 2tanh^^(\/^), and is depicted as a black curve 
in Fig. 21 It provides an upper bound for all two-mode 
Gaussian states distributed in the {E'j^,A'^} plane, in 
the region of moderate entanglement; such a bound con- 
verges to 1 in the separability limit {iy — 1). We can thus 
conclude that the degree of nonclassicality of correlations 
in separable Gaussian states is always nonzero (apart 
from the trivial case of uncorrelated, product states) but 
stays nevertheless limited: It can at most reach unity, 
whether measured by the discord [13, [11] or by the Gaus- 
sian AMID. 

The second expression in Eq. (j66p corresponds instead 
to the Gaussian AMID of the states of Eq. (p^ in the 
limit a — > oo. It bounds from above the value of A'~^ 
for all two-mode Gaussian states with fixed i?^ ^ 0.441 
[where this number is obtained by setting i> — 4/e — 1 in 
Eq. In the limit of infinite entanglement (i> — )■ 0), 

the upper bound on the Gaussian AMID converges to 
E"^ + In 4 — 1 . Combining this observation with the lower 
bound (pT|) . we have that, interestingly, the following 
sandwich relation holds for all two-mode Gaussian states 
with Ef :^0, 

EfipAB) < A'^ipAB) < EfipAB) + ln4 - 1 . (67) 

In the previous analysis, we have identified some sim- 
ilarities as well as some key differences in the quantifica- 
tion of nonclassical correlations versus entanglement of 
Gaussian states when Gaussian AMID rather than quan- 
tum discord are employed. In order to have a visual 
comparison between the two nonclassicality indicators 
and the Gaussian EoF, we focus on the relevant two- 
parameter class of symmetric squeezed thermal states 
Pab' with CM as in Eq. (15^ . For these states, the entan- 
glement is given by Eq. (j64p (independently of a), while 
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states, E'j can be exactly recast as a function of T)*^ 
and JP: Knowledge of two nonclassicality quantifiers 
fixes the third one univocally. We remark that such a 
strict result does not extend to more general two-mode 
Gaussian states, which distribute filling a more com- 
plex, finite- volume three-dimensional region in the space 
{T>^,A^, Ef}. Figures [I(c)] and|lof this paper, and Fig- 
ure 1 (right) of Ref. yjj], represent the two-dimensional 
projections of such a region onto the planes {V^jA^}, 
{Ef,A^}, and {Ef,!)^}, respectively. 



VII. CONCLUSIONS 



FIG. 5: (Color online) 3D Plot of Gaussian EoF Ef [Eq. ^ ] 
versus Gaussian discord [Eq. (|68p ] and Gaussian AMID 
[Eq. (|69|) ] for two-mode symmetric squeezed thermal 
Gaussian states, characterized by their covariance parame- 
ters a and D [see Eq. (|62|) ] . The solid black line accommodates 



pure states {i> - 
dimensionless. 



a — y/a'^ — 1) . All the quantities plotted are 



the discord can be written as llTl \li 



2(1 + a) 



(4a(i> + 1) - 2i>^) tanh"^ 



2aD 



4(a + l)v/^(2a^^tanh ^ 



In 



a- 1 



In 




(a-l)(i>+l)(2a-:>+l) 

(68) 



and the Gaussian AMID reads 
A^ipTB) = -hi 



2av — — 1 



a? - 1 



2acoth ^(a) 



2^Jv{2a~v) tanh"^ 



1 



I (a+l)^-(a-j>)^ ' 



^ i>(2a — v) j 
1 + a (4 + a (4 - 2ai> 

otherwise. 



^))>0; 



(69) 



Fig. [5] shows Ef plotted versus I?^ and A'^ for this 
particular class of Gaussian states, spanned by a and 
V. All two-mode symmetric squeezed thermal states sit 
on a two-dimensional surface in the space of the three 
entropic nonclassicality indicators, providing a direct ev- 
idence of the intimate yet intricate relationship between 
the different aspects of quantumness in Gaussian states. 
One can notice the branch of separable states, in the 
plane Ef = with generally nonzero discord and Gaus- 
sian AMID; while for all entangled squeezed thermal 



We have performed an exhaustive study of nonclas- 
sical correlations in generic two-mode Gaussian states 
using information-theoretic nonclassicality quantifiers, in 
particular the 'measurement-induced disturbance' (MID) 
and its Gaussian-optimized version, the Gaussian 'ame- 
liorated measurement-induced disturbance' (Gaussian 
AMID). For a given Gaussian state the MID is a gap 
between its quantum mutual information — quantifying 
the total correlations — and the classical mutual infor- 
mation of outcomes of local Fock-state detections (detec- 
tions of local eigenprojectors) — that captures a specific 
type of non-Gaussian classical correlations in the state. 
The Gaussian AMID is, on the other hand, a gap between 
the quantum mutual information and the maximal clas- 
sical mutual information that can be obtained by local 
Gaussian measurements, the latter quantifying the max- 
imum classical correlations that can be extracted from 
the state by local Gaussian processing. An analytical 
form of the Gaussian AMID can be derived for the impor- 
tant subclasses of symmetric states and squeezed thermal 
states which include pure states, while for a generic mixed 
Gaussian state one has to find roots of a higher-order 
polynomial in a single variable, which can be solved effi- 
ciently by numerical means. Further analysis reveals that 
MID is mostly larger than Gaussian AMID and therefore 
overestimates the amount of nonclassical correlations. In 
fact, for a fixed value of Gaussian AMID, it is possi- 
ble to find states with an arbitrarily large MID, even 
if the Gaussian AMID is infinitesimally small. On the 
other hand, there also exists a volume of Gaussian states 
encompassing pure states for which the MID is strictly 
smaller than the Gaussian AMID, which surprisingly un- 
veils the importance of non-Gaussian measurements for 
the correct assessment of the amount of nonclassical cor- 
relations in Gaussian states. 

We have further compared the MID and Gaussian 
AMID with the two-way Gaussian discord. We found 
that again there is no upper bound on MID for a fixed 
value of discord but a close upper bound (as well as a 
lower one) does exist for the Gaussian AMID. Finally, 
we have also compared the Gaussian AMID with the 
Gaussian entanglement of formation, identifying lower 
and upper bounds for the former as a function of the 
latter. In particular, the Gaussian AMID turns out to 
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always exceed the Gaussian entanglement of formation 
for all two-mode Gaussian states, enforcing a novel hier- 
archy between two different forms of nonclassicality. Ex- 
act relations between Gaussian AMID, Gaussian discord 
and Gaussian entanglement of formation can be formu- 
lated for special families of Gaussian states, such as the 
symmetric squeezed thermal states. 

On a more technical side, we have shown that symmet- 
ric measures of (non)classicality of correlations such as 
the classical mutual information and the Gaussian AMID 
are tractable by analytical tools for important subclasses 
of states even considering optimization over local gener- 
alized Gaussian measurements. Our results also demon- 
strate that non-Gaussian processing for correct quantifi- 
cation of (non)classical correlations in Gaussian states is 
in order. 

We believe that these findings will inspire further re- 
search on the characterization of quantum correlations in 
the Gaussian scenario and beyond. 
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Appendix A: Reduction to covariant rank-one 
POVMs 

Here we prove that covariant rank-one POVMs, i.e. 
POVMs of the form Eq. ([5]) with a pure seed state Ilj, are 
optimal among all Gaussian POVMs for the evaluation of 
the classical mutual information or, equivalently, of the 
Gaussian AMID, Eq. §EI- 

The measurement of local POVMs ([5]) on a Gaussian 
state pAB gives the outcome d — {d^,d]^)^ distributed 
according to 

P{d) = TT[IlA{dA) ® flB{dB)pAB]- (Al) 

Making use of the overlap formula for Wigner functions 
(isj , the distribution can be expressed as 

P{d) = {27r)' J W^^^,^^{rA)W^^^,^^{rB) 

xWp^^{rA,rB)d^rAd\B- (A2) 

Substituting into the formula from Eq. ^ , where we set 
r = (r^i, r^y^ , and performing the integration, we obtain 



the distribution (jAl[) in the form: 

P{d) = ^ ^-d-^ (1+1 A®iB)-'d ^^3) 

7r2^det(7 + ^a®1b) 

where 7_a.s are CMs of the seed elements of POVMs ^ 
and 7 is the CM of the state pab- The CMs 7_a.s can 
be expressed as 7j = ^^^^ where j^^^ — S^^{Sj)^^ 

is a pure-state CM {Sj symplectically diagonalizes jj) 
and Nj = (vj — 1)7^'^'' is a positive-semidefinite matrix 
{vj > 1 is a symplectic eigenvalue of 7^). Therefore, the 
outcome dj, j = A, B oi a. generic POVM Ilj{dj), with 
the seed element Ilj being a mixed state with CM jj, 

can be expressed as dj = d'^'^ + Xj, where dj'^^ is the 
outcome of POVM with pure-state seed element having 
CM 7j and xa^Xb are mutually uncorrelated random 

variables uncorrelated with d^^d^'' obeying Gaussian 
distributions with classical correlation matrices Na and 

(tt) 

Nb, respectively. Since such processing of variables dj 
cannot increase their Shannon mutual information due to 
what is now known as data processing inequality (derived 
first for continuous random variables in [46|), we can re- 
strict without loss of generality to optimization over pro- 
jections onto pure states. 

Note, that covariant measurements ([5]) with pure- 
state seed elements maximize classical mutual informa- 
tion even within the framework of a larger class of gener- 
ally noncovariant Gaussian POVMs possessing the struc- 
ture [H, 113 

^.(^.) = ^^D,(d,)fl,(y,)D]{d,), J = A, B, (A4) 

and satisfying the completeness condition Ttj{zj)dzj — 
1 j . Here pj (yj ) is a normalized distribution of the param- 
eter yj, Ilj{yj) is a normalized Gaussian state with CM 
"ij{yj) dependent on parameter yj and Zj — {dj,yj)'^. 
Upon measuring the POVM ljA4[) on the Gaussian state 
Pab 1 one finds the outcomes za and zb to follow the dis- 
tribution VizA.ZB) = pA{yA)PB{yB)P{d,yA,yB), where 
the distribution P(d,yA,yB) is obtained from Eq. (|A3p 
by replacing 7^ with ^jijjj)- Denoting the classical 
mutual information of the distribution 'P{zatZb) and 
Pid,yA,yB) as I[A{za) : B{zb)] and I[A{yA) : S(ys)], 
respectively, one then has 

I[AizA) : B{zb)] = f I[A{yA) : B{yB)] H p{yj)dyj. 

J j=A,B 

(A5) 

Hence it follows immediately that 

IciPAB) < , max I[AizA) : BizB)] 

= T[A{y\) : B{y%)], (A6) 

where we have to maximize over all Gaussian POVM 
elements T^j{zj) and y^, j = A,B label the POVM ele- 
ments Hjiy'j), which maximize T[A(jja) ■ B(jjb)]- If we 
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take these as seed elements of POVMs ([5]), we construct X[A{y'^) : B{y'^)] and therefore achieve the classical mu- 
local covariant POVMs which give mutual information tual information I^(/5ytB)- 
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